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Dark-matter sterile neutrinos in models with a gauge singlet in the Higgs sector 
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Sterile neutrino with mass of several keV can be the cosmological dark matter, can explain the 
observed velocities of pulsars, and can play an important role in the formation of the first stars. 
We describe the production of sterile neutrinos in a model with an extended Higgs sector, in which 
the Majorana mass term is generated by the vacuum expectation value of a gauge-singlet Higgs 
boson. In this model the relic abundance of sterile neutrinos does not necessarily depend on their 
mixing angles, the free-streaming length can be much smaller than in the case of warm dark matter 
produced by neutrino oscillations, and, therefore, some of the previously quoted bounds do not 
apply. The presence of the gauge singlet in the Higgs sector has important implications for the 
electroweak phase transition, baryogenesis, and the upcoming experiments at the Large Hadron 
Collider and a Linear Collider. 

PACS numbers: 14.60.St, 95.35.+d 



I. INTRODUCTION 

The discovery of the neutrino masses can easily be 
incorporated into the Standard Model (SM) by adding 
two or more SU(3)xSU(2)xU(l) singlet fermions, often 
called right-handed neutrinos, which are allowed to have 
the Yukawa couplings to the Higgs boson and the stan- 
dard, left-handed neutrinos. The Yukawa couplings gen- 
erate the Dirac mass terms for the neutrinos after the 
spontaneous symmetry breaking. In addition, the singlet 
fermions can have some Majorana masses. The interplay 
between the Dirac mass and the Majorana mass, known 
as the seesaw mechanism [l|, can accommodate the ob- 
served neutrino masses for a variety of Majorana masses. 
If the Majorana mass terms are large, the particles asso- 
ciated with the singlet fields are very heavy. However, if 
one or more Majorana masses are below the electroweak 
scale, the so called sterile neutrinos appear among the 
low-energy degrees of freedom. These new particles can 
be the cosmological dark matter 0, d, 13, S 0, S Hi , their 
production in a supernova can explain the pulsar kicks [oj 
and can affect the supernova explosion in a variety of 
ways 10:]; the same particles can play an important role 
in the formation of the first stars and some other 
astrophysical phenomena [T^ . 

The properties of the sterile dark matter, and, in par- 
ticular, how warm or cold it is for a given mass, de- 
pend on the production mechanism. One mechanism, 
which generates a population of relic sterile neutrinos 
at the sub-GeV temperature was proposed by Dodel- 
son and Widrow (DW) [2]. If the lepton asymmetry is 
negligible, this scenario appears to be in conflict with a 
combination of the X-ray bounds and the Lyman-a 
bounds [13, [l^ . This conclusion is based on the state- 
of-the-art calculations of the sterile neutrino production 
in neutrino oscillations It is possible to evade this 
constraint if the lepton asymmetry of the universe is 
greater than 0(10"'^) \M- O n the other hand, some as- 
tronomical observations |16l . Iv^ point to a non-negligible 
free-streaming length for dark matter, which favors warm 
dark matter. Moreover, warm dark matter can cause fil- 



amentary structure on small scales |18l | , in contrast with 
cold dark matter. It is also possible that the sterile neu- 
trinos make up only a fraction of dark matter [l^, in 
which case they can still be responsible for the observed 
velocities of pulsars 0, 0] ■ 

Dodelson-Widrow mechanism is not the only mecha- 
nism by which sterile dark matter could be produced. 
The relic population of sterile neutrinos could be gener- 
ated in a variety of ways, for example, from their cou- 
pling to the inflaton Q, the electroweak-singlet Higgs 
boson [tI, or the radion Whatever the production 
history of sterile neutrinos might be at the high tem- 
perature, there is always some additional amount pro- 
duced in neutrino oscillations at some sub-GeV temper- 
atures [3, The two components can have very dif- 
ferent momentum distributions. Therefore, generically 
this form of dark matter is a mixed two-component dark 
matter, which can have some very interesting observable 
consequences jl9| . 

In this paper we concentrate on the possibility that the 
relatively light Majorana mass could arise via the Higgs 
mechanism in a model with an SU(2) xU(l)-singlet Higgs 
boson coupled to the Standard Model Higgs boson [20| . 
and that the sterile neutrinos could be produced from the 
Higgs decays at a temperature as high as 100 GeV Q- We 
will explore various scenarios for such production and the 
implications for the electroweak phase transition. In par- 
ticular, we will address the cooling and the red-shifting of 
dark matter, which have important implications for dark 
matter profiles in halos [l^ [l3| , the small-scale structure 



inferred from Lyman-a observations [1^ 
ity dispersion in dwarf spheroids [T6| . 



and the veloc- 



II. MAJORANA MASSES FROM AN 
EXTENDED HIGGS SECTOR 

Although the Standard Model was originally formu- 
lated with massless neutrinos i^i transforming as compo- 
nents of the electroweak SU(2) doublets {a = 1, 2, 3), 
the neutrino masses can be accommodated by a relatively 
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minor modification. One adds several electroweak sin- 
glets Na (a = 1, ...,n) to the Standard Model and builds 
a seesaw lagrangian [Ij: 

£ = CsM+lNama~yc.aH^ Lo^N^-^^ N^Na + h.C. (1) 

The neutrino mass eigenstates (i = 1, ...,n + 3) 
are linear combinations of the weak eigenstates {va, Na}. 
They are obtained by diagonalizing the mass matrix: 



yo.a{H) 
yaa{H) diag{Mi,...,M„} 



(2) 



As long as all yaa{H) <C Ma, the eigenvalues of this ma- 
trix split into two groups: the lighter states with masses 
of the order of y'l^{H)'^ /Ma, and the heavier eigenstates 
with masses of the order of Ma- As usual, we will call 
the former active neutrinos and the latter sterile neutri- 
nos. The mixing angles in this case are of the order of 

The number n of the right-handed singlets is unknown, 
although it is clear that n > 2 is a necessary condition to 
explain the results from the atmospheric and solar neu- 
trino experiments [2T| . Theoretical considerations do not 
constrain the number n of sterile neutrinos. In particular, 
there is no constraint based on the anomaly cancellation 
because the sterile fermions do not couple to the gauge 
fields. The experimental limits exist only for the larger 
mixing angles [1^. The scale of the right-handed Majo- 
rana masses. Ma, can vary over many orders of mag- 
nitude. It can be much greater than the electroweak 
scale or it may be as low as a few eV [1^. It is 
also possible that some of the right-handed Majorana 
masses are much larger than others. The seesaw mech- 
anism can explain the smallness of the neutrino masses 
even if the Yukawa couplings are of order one, as long as 
the Majorana masses Ma are large enough. However, the 
origin of the Yukawa couplings remains unknown. If the 
Yukawa couplings arise as some topological intersection 
numbers in string theory, they are generally expected to 
be of order one 12711 . although very small couplings can 
are also possible [2g. However, if the Yukawa couplings 
arise from the overlap of the wavefunctions of fermions 
located on different branes in extra dimensions, they can 
be exponentially suppressed and are expected to be very 
small [2^. If one or more singlets have Majorana masses 
below the electroweak scale, they can appear as sterile 
neutrinos and can have important ramifications; for ex- 
ample, dark matter can be made up of sterile neutrinos 
with mass of several keV [2], and the same particle can 
be responsible for the observed pulsar kicks [9| . 

Several recent papers have studied in detail one partic- 
ular case, named i^MSM [5], which corresponds to n = 3, 
Ml ~ keV, and M2 « M3 ~ 1 - 10 GeV. In this model, 
the keV sterile neutrino serves as the dark matter particle 
(and can explain the pulsar kicks), while the degenerate 
heavier states, M2 ~ M3, make the model amenable to 
leptogenesis by neutrino oscillations 30] . 



The possible role of keV sterile neutrinos in astro- 
physics and cosmology, from dark matter to pulsar kicks, 
to early star formation, makes the possibility of their exis- 
tence very intriguing. However, if the neutrino Majorana 
masses Ma are below the electroweak scale, one should 
try to explain the origin of this scale. The other fermions 
in the same mass range acquire their masses from the 
Higgs mechanism. Can the mass terms in eq. ([1]) also 
arise from the Higgs mechanism? The answer is yes; this 
requires an extension of the Higgs sector by an SU(2) 
singlet field coupled to the righted-handed fermions as in 
Rcfs. fi,0,[2|: 



C = CsM + iNa^Na - yc^aH^ L^Na 
- V{H,S) + h.c. 



fa 



S KNa 



(3) 



We will assume that S' is a real scalar field to avoid 
the light Nambu-Goldstone bosons associated with the 
breaking of the lepton number U(l); the presence of such 
light bosons would render the sterile neutrinos unstable, 
hence they could not be dark matter (although they could 
still explain the pulsar kicks [9]). If the singlet has very 
small mass and a large VEV, it can be the inflaton [y]. 
We will not discuss this interesting possibility here, but 
we will concentrate instead on a singlet Higgs whose mass 
and VEV are both of the order of 100 GeV, which, inci- 
dentally, is the requirement for the keV dark matter, as 
long as the mass and VEV of S are of the same order of 
magnitude 7]. 

As soon as the SNN coupling is introduced in the la- 
grangian, there appears a new way in which the relic pop- 
ulation of sterile neutrinos can be produced, namely from 
the decays S — > NN . This decay mechanism can oper- 
ate in addition to the neutrino oscillations mechanism of 
Dodelson and Widrow and one has to compare the 
relative amounts produced by each of them. Another 
important issue is how cold the dark matter is if it is 
produced predominantly from the Higgs decays. Since 
the production occurs mainly at temperatures of the or- 
der of the Higgs mass, T ~ 100 GeV, the reduction in 
the number of degrees of freedom and the entropy pro- 
duction that takes place as the universe cools down from 
T ^ 100 GeV causes the dark matter population to be di- 
luted and red shifted by a factor ^ > 33 in the density and 
factor > 3.2 in the average momentum. These values 
reflect only the Standard Model degrees of freedom, and 
any additional new physics will make ^ even larger. The 
corresponding free-streaming length is shorter, and the 
Lyman-alpha bounds become proportionately weaker [3] . 

In this paper we discuss the details of sterile dark mat- 
ter production in a model represented by the lagrangian 
([3]), with the scalar potential 

V(H,S) = -/4|Hp-i452^ia53+c.|i7p5 

+ \h\H\'' + ^\sS^ + 2Xhs\H\''S'' (4) 
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III. STERILE DARK MATTER: COLD OR 
WARM? 

If dark matter has a non-zero free-streaming length, 
the structure on small scales may be suppressed. Stud- 
ies of small-scale structure based on the observations of 
dwarf Spheroids [l6i| or Lyman-a forest data Ti*] can con- 
strain or measure the free-streaming length of the dark- 
matter particles, but the relation between this length and 
the particle mass depends on the production mechanism. 
One can approximately relate the free-streaming length 
to the mass rrig and the average momentum of the sterile 
neutrino: 

This is a relatively good measure in many cases, although 
in general one has to calculate the full power spectrum. 
The observations of Lyman-a forest constrain the free- 
streaming length to be less than 0.11 Mpc This 
bound does not translate directly into a constraint on 
the mass because the average momentum depends on the 
production mechanism. For three scenarios usually dis- 
cussed in the literature, 

/ / \ \ ( 0.8-0.9, forDW 

( ^ j = w 0.6, for L 7^ 0, resonance 

V^-^^^^T«koV [ <0.2, for Tprod > 100 GeV 

(6) 

Here DW stands for Dodelson-Widrow production mech- 
anism via non-resonant neutrino oscillations Q , 7^ 0" 
refers to the Shi-Fuller production via the resonant neu- 
trino oscillations in the case when the lepton asymmetry 
is relatively large [4(], and "Tprod ^ 100 GeV" refers to 
the production of sterile neutrinos at a temperature well 
above the QCD scale, in which case the cooling and re- 
duction of the degrees of freedom causes the red shift in 
the population of dark matter [3] . 

For the same mass, the sterile dark matter can be 
colder or warmer, depending on the production mecha- 
nism. This is clear from equations ([5]) and ([6]), which, for 
a given cosmological scenario, relate the free-streaming 
length with the mass. Therefore, we will pay close atten- 
tion to the factors that can affect the momentum distri- 
bution in each scenario. 

There are several ways in which the population of dark 
matter particles could have formed in our model: 

• The bulk of sterile neutrinos could be produced 
from neutrino oscillations. If the lepton asymmetry 
is negligible, this scenario appears to be in con- 
flict with a combination of the X-ray bounds [T3| 
and the Lyman-a bounds [T5| . although it is pos- 
sible to evade this constraint if the lepton asym- 
metry of the universe is greater than O(10~^) [3]. 
It is possible that the decays of additional, heavier 
sterile neutrinos, can introduce some additional en- 
tropy and contribute to cooling of dark matter . 



It is also possible that the sterile neutrinos make up 
only a fraction of dark matter 0, [lB| , in which case 
they can still be responsible for the observed veloc- 
ities of pulsars. 

• The bulk of sterile neutrinos could be produced 
from decays of S bosons at temperatures of the 
order of the S boson mass, T ~ 100 GeV. This sce- 
nario was discussed in Ref. 0]. In this case, the 
Lyman-a bounds on the sterile neutrino mass are 
considerably weaker than in the former case. 

• The decays described above could happen before 
a first-order phase transition, and the entropy re- 
lease in the transition could redshift the popula- 
tion of the dark-matter particles. We have explored 
this possibility in detail, as discussed below, but we 
have not found a range of parameters in which the 
phase transition could cool down the sterile dark 
matter significantly. 

• S bosons could be so weakly coupled to the rest of 
the Higgs sector that they would go out of equi- 
librium and decay out of equilibrium at some tem- 
perature T < 100 GeV. As discussed below, this 
scenario can produce a sufficient amount of dark 
matter. 

We will now discuss these possibilities in detail. 

IV. PRODUCTION FROM THE HIGGS 
DECAYS IN EQUILIBRIUM 

The interactions of the singlet Higgs bosons with SM 
particles have been studied by McDonald in Ref. [11], 
where the S bosons were made stable by imposing a 
global U(l) symmetry, which removed the odd power cou- 
plings, and by setting < 0, which forced {S) = 0. In 
this case, the coupling Xhs controls the SS XX an- 
nihilations, into SM fermions and the W, Z bosons. We 
do not require S to be stable. After S develops a VEV, 
other couplings also contribute to the annihilations into 
SM particles. For each of these processes the cross sec- 
tion for annihilation is: 

aann-^10-^% (7) 

At some temperature, these processes fail to keep the 
S particles in equilibrium, and they freeze out at Tf = 
ms/rf. For very small Xhs ^ 10~^, S bosons never come 
into equilibrium. A more detailed numerical calculation 
yields the dependence of the freeze-out time parameter 
Tf on A^f5 shown in fig. [1] 

The cubic couplings contribute to the annihilation pro- 
cesses through exchange of virtual S bosons. In fact, this 
will be the dominant process that keeps S particles in 
equilibrium, as long as ^ > Xhs, where ms is the S 
boson mass. Comparing with Fig. [H one can see that 
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FIG. 1: The variation of the S boson freeze-out parameter 
n = ms/Tf with the coupling to SM particles \hs- For 
numerical estimations we used ms = 200 GeV. 



this process can keep S bosons in equilibrium down to 
rather low temperatures, even if a,Lu are well below the 
S mass. 

Let us now assume that the a, lu and/or Xhs are large 
enough (exact limit for Xhs to be defined below) to keep 
S in equilibrium down to temperatures well below its 
mass. One can make a rough estimate of the sterile neu- 
trino production by multiplying the S number density 
by the S NN decay rate F = f^ms / {16n) and by the 
time available for the decay, r ~ Mo/2T^, at the latest 
temperature at which the thermal population of S is still 
significant, namely T ~ ms- At lower temperatures, the 
S number density is too small, much smaller than T^. 
One obtains an approximate result 



T3 



r 



Mo 

2^2 



IGtt ms ' 



(8) 



1/2 



10^** GeV is the reduced 



where Mo = (^f^^ 
Planck Mass. 

This simple estimate is in agreement with the solution 
of the kinetic equation discussed below. Of course, this 
description breaks down if the S particles decouple and 
decay at a much later time. We will come back to this 
possibility. For now, let us assume that S particles main- 
tain their equilibrium populations down to temperatures 
at least a factor of a few below their masses. The dark 
matter abundance for sterile neutrinos from the decay of 
bosonic particles in equilibrium was first computed for a 
model in which the S field served as the inflaton with a 
potential adjusted to have (S) ^ ms The results of 
this computation carry over to our case. Here we do not 
require S to be the inflaton, and we take (S) ^ ms, as in 
Ref. [rj. As was shown in Ref. [t'I, the choice {S) ^ ms, 
along with the requirement that sterile neutrinos make up 
all the dark matter, force the S boson mass and VEV to 
be right at the electroweak scale, (S) ~ ms ~ 10^ GeV, 
suggesting that S may, indeed, be a part of the Higgs 
sector of the extended Standard Model, and justifying 
some of our assumptions regarding the Higgs potential. 



Dark matter production from particle decays has been 
considered in a number of papers [13]. Let us first 
consider the decaying particle in equilibrium. As in the 
case of the inflaton decay ^ , the sterile neutrino distri- 
bution function n(p, t) is found from the following kinetic 
equation: 



dn ^ dn 
dt dp 



2msT 

p2 



nsdE, 



(9) 



where F — msp /IQtt is the partial width of the S boson 
decay. It is assumed that the sterile neutrinos are never 
in equilibrium, and the inverse decays NN — > S can be 
neglected, which is true for small Yukawa couplings / < 
10~^. Transforming to the variables: r = x 
one can rewrite eq. ([9]): 



p. 

T' 



dn _ p Mo r2 
dr Stt 



X + -, 



(10) 



Since S and H mix, one has to consider the mixed mass 
eigenstates in plasma. Both of them can decay into sterile 
neutrinos. The SM Higgs is in thermal equilibrium due to 
the coupling with the SM particles. For the temperature 
range in which S is also in equilibrium, the distribution 
functions of the two mass eigenstates will be: 



1 



(11) 



Then eq. ^ yields 



n^{x, r) 



y^pMo 
^ Stt to, 



+ 



l+lfr (z-l)3/2dz 



1 



(12) 



where the subscript j = 1,2 runs over the two Higgs 
mass eigenstates and the superscript O denotes produc- 
tion from decays of S bosons in equilibrium. (We will use 
for the case of S bosons decaying out of equilibrium.) 
In (|12p the first term is important when r < 1, while the 
second is the dominant one for r > 1. The above solu- 
tion was obtained assuming fj, mj and the number of 
degrees of freedom (7* remain constant. This is not valid 
after the electroweak phase transition takes place, since 
the Higgs mass eigenvalues and their mixing are different 
in the two vacua. If fj , mj or change at some points 
in the evolution of the universe, the solution has to be 
adjusted to include the contributions from all the periods 
corresponding to different fj , mj , . Each of these con- 
tributions will still be given by (|12p . for the appropriate 
values of the parameters and taken over the respective 
time intervals. 

This complication turns out to be irrelevant, since the 
production rate of sterile neutrinos through each mode 



mj exhibits a peak at rj 
duction temperature 



2.3, which defines the pro- 
: mj/2.3. Most of the sterile 
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neutrinos are produced around that temperature, and by 
the time when rj « 10 the production of sterile neutrinos 
though rrij decays has practically been completed. More 
specifically, in the simplified case mi ^ TO2, /i 3> /2, 
TOi decays will dominate over m2 decays and the total 
abundance of sterile neutrinos K, = Ns/ s at any later 
temperature will be: 



45 



Mn 



327r5 g*(Tpi.od) m 



Vir) 



(13) 



where 



2/W 



da; 



4n(l-e"^"5:^) 



(Z- l)3/2dz 



1 



(14) 



and we have dropped the mass eigenstate index for sim- 
plicity. The 5*(rpi.od)~"'^ factor in (fT5)) designates the fact 
that the sterile neutrino population will be diluted by 



C = 5*(rprod)/5*(0.1MeV) 



(15) 



as the universe cools down, due to the entropy release as 
the effective degrees of freedom decrease. At r ^ oo, the 
sterile neutrino abundance produced from in-equilibrium 
decays takes the limiting value: 



27C(5) p 



Mn 



(16) 



We require that the decays of S bosons occur in equi- 
librium, that is Tf < j^. Then we get from Fig. [1] 
^HS ^ 5 • 10"^, for a = w = 0. However, if a, u; > 1 GeV 
S bosons stay in equilibrium down to the desired tem- 
perature, regardless of Xhs- 

The momentum distribution in eq. (|12p is non-thermal. 
Taking into account only the dominant decay mode, one 
obtains (same as in Ref. 6]) the momentum distribution 
of dark-matter particles at r — > oo: 



Stt to 



(Z- l)3/2dz 



1 



(17) 



for which the average momentum at temperature 
T ^ 10^ GeV, immediately after their production, 
is (c/. Ref. i) 



M 

T 



T~WO GeV 



378 C(5) 



2.45. 



(18) 



This is lower than the same quantity for a thermal dis- 
tribution, {p)/T = 3.15. 

Even more importantly, these momenta are further 
redshifted as the universe cools down from the temper- 
ature at which most dark matter is produced, Tprod ~ 
TO 100 GeV to the much lower temperatures at which 
the structure begins to form. As the universe cools down. 



the number of effective degrees of freedom decreases from 
g*{Tprod) = 110.5 to 3^0.1 MeV) = 3.36. This assumes 
no new physics below the Higgs mass; any new physics 
would cause an additional cooling of the dark matter. 
The ratio of dark matter to entropy is reduced by the 
factor ^ « 33. This causes the redshifting of (ps) by the 
factor ^1/3: 



0.76 



110.5 



1/3 



(19) 



/ (T«;lMoV) 

This is very different from the DW scenario , in which 
the average neutrino momentum at low temperature T is 



'DW 



2.83 T. 



(20) 



Comparing eqns. ([^0]) and ([TO]) , one concludes that the 
sterile neutrino mass corresponding to the same free- 
streaming length can vary by more than a factor of 3 
depending on the production scenario A detailed 
analysis of the free-streamingproperties of "chilled" dark 



matter is presented in Ref. [34 1 



The dark matter abundance in this model depends 
on the details of the Higgs mass matrix and the two- 
component decays, aside from which it has the form: 



fl^, 0.2 



/ 



10- 



{S) 



33 



(21) 



Since we expect the masses of the two mass eigenstates to 
be of the same order, and considering the cubic power of 
the unknown coupling /, the details of the solution are 
not very important. However, what may be important 
is the additional effect of the first-order phase transition 
on the average momentum of the dark matter particles. 
If the dark matter population is redshifted significantly 
by the entropy release in the phase transition, then the 
Lyman-a bound could be further relaxed. 



V. ELECTROWEAK PHASE TRANSITION 

As discussed above, the population of sterile neutri- 
nos is subject to dilution and redshift due to the entropy 
production that occurs (i) in any possible phase transi- 
tions at lower temperatures, and (ii) when the number of 
degrees of freedom in plasma decreases due to the decou- 
pling of Standard Model particles below 100 GeV. The 
redshift due to (ii) alone can reduce the momenta of dark 
matter particles by a factor more than 3.2 Of course, 
the dilution due to (i) matters only if a first-order phase 
transition takes place after the sterile neutrinos are pro- 
duced, and, moreover, if S bosons are too heavy to have 
a high number density in the new vacuum after the phase 
transition. To study this possibility, one has to take into 
account the temperature effects on the effective poten- 
tial and the history of the phase transitions predicted by 
the model. Electroweak phase transition in a model with 
a singlet Higgs has been analyzed in Refs. [1^ [s^, [IJ . 
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The plausibility of the first-order phase transition makes 
the electroweak baryogenesis a viable possibility, and it 
has implications for the LHC and the ILC [s^]- Here we 
concentrate on the effects the first-order transition could 
have on the population of dark-matter sterile neutrinos. 



The temperature-dependent contribution to the effec- 
tive potential at one loop is [H, [s^ : 



rp2 



(27) 



A. Finite-temperature effects and the first-order 
piiase transition 



where 



hj{y) 



dxx log 



(28) 



The tree level effective potential in terms of the VEV 
of the two Higgs bosons, (H) — -^rj and (5) = ct is 



^ 2 2 



A 



^ ,,2 ^2 , "^s 4 



(22) 



To study the phase transition, we included the one-loop 
temperature-dependent corrections and analyzed the po- 
tential numerically, as discussed below. 

The tree level Higgs mass eigenvalues are: 

KJ' = i[(3A,+2A„,)7,2^(3A,+2A„Ja2 

+ {u; + a)a - fil - nl 

± {[{3X, - 2X„,)^^ - (SX, - 2X,,)<j' 

+ (w-a)cr-/4 ^_^2j2 

1/2' 



+ Arj^AX,,a + ujf} 



(23) 



The 1-loop, zero temperature correction to the effective 
potential [315 in the MS renormalization scheme is: 



The upper sign corresponds to bosons, while the lower 
one to fermions. Since we consider a wide range of tem- 
peratures, we do not make use of the well known high-T 
expansion of the functions Vfhig ^hc ring contri- 

bution of the gauge, Higgs and goldstone bosons: 



T 
12 



— |Tr (to 

TT L L 



H„ 



^3/2 



Tr 



gb T ^^cjb) - yi'I'gb 

TT ^,3/2 _ f 2 \3/2 
-l^-lhiggs; v"-higgs/ 



(m 



H,)3/^ ~ (m^)3/2 



(29) 



where mhiggs is the tree level Higgs mass mixing matrix, 
corresponding to the potential ([22|) . whose eigenvalues 
are given in (|23p. TOgj^ is the mass mixing matrix for the 
electroweak gauge bosons: 



(30) 



/ 


2 2 

9 V 
4 








\ 







4 


















gg'v^ 




4, 2 


,2*2 


I 








gg V 

4 


g V / 

4 / 



Hi are the thermal contributions to the masses, given for 
our model byHHHIJ: 



T.i/ X 4/ N /', m,(?7,cr)2 3 



(24) 

where Ui are the degrees of freedom of the contributing 
particles and mi{r],a) are their field-dependent masses. 
The main contributions are from the neutral component 
of the SM Higgs, the singlet Higgs, the Goldstone bosons 
X, the W and Z gauge bosons and the top quark t: 

n-t = -12, 7i„, = 6, = 3, = 3, n„ = = 1 

(25) 

The gauge-singlet Higgs S does not couple to the 
fermions or the gauge bosons, thus their field dependent 
masses are the same as in the minimal SM : 



m 



^rf, to2 
2 ' 



9 2 2 
77 , TO 

4 ' ' z 



9' 



- o'2 



(26) 



The Higgs mass eigenvalues are given by mi{ri^, ctq) 
stand for particle masses at the vacuum state (770, do) at 
zero temperature. 



H 



gb 



Hhiggs — 



diag 
diag 



11 



g^T\ ^-9'T\ ^-g'T\ ^-g'-'T 



11 



11 



11 



,/2rp2 



16^ 16^ 

1, 4 
4A, + -A„, 



16- 



1 

16 



A, 



9' 



yt_ 
4 



4 



3 



'HS 

3 



(31) 



The effective potential at finite temperature is the sum 
of (El, (El and (EZl): 

K=ff(?7,'T,T) =F°(?7,a)+Fi(77,a)+F^(77,(T,T) (32) 

In the calculations that follow we ignore for simplicity, 
since it is only a small correction to the zero-T effective 
potential, and we treat the imaginary part of the poten- 
tial as usual [ilj . 

The history of the universe for a typical set of param- 
eters discussed in table HI is shown in fig. El Due to the 
symmetry 77 — > —77, only the 77 > half plane need to be 
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T=Tc 
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FIG. 2: Contour plots of V^ftiv, o". T) at T > Tc, T = Tc, and 
T — 0, corresponding to the parameter set A of table HI At 
T ^ Tc the universe is in the unbroken phase rj — 0, with 
o" > 0. At lower temperatures, the minimum of Kft shifts to 
non-zero i], while a second minimum appears in the a < 
region. The two minima become degenerate at T = Tc, and 
at T = the true vacuum is located at r; = 246 GeV, a < 0. 



considered. Since the singlet does not couple to the gauge 
bosons and the t quark, it receives the smallest correc- 
tion at a high temperature. Therefore, there is usually a 
range of temperatures in which the Higgs doublet has no 



VEV, while the singlet has a VEV. At a higher temper- 
ature, the singlet VEV also tends to 0, although it never 
completely disappears due to the ct — s- — cr asymmetry, 
induced by the a and oj terms. 

Therefore, in the early universe, at T ^ 100 GeV, the 
effective potential has a unique minimum at {0,crf}. At 
a lower temperature, for some range of the parameter 
space, the doublet also develops a VEV and the universe 
shifts to {rif{T),af{T)} through a second order phase 
transition. In the meanwhile, a local minimum has been 
developed at {?7t(T), crt(r)}. At T = Tc the two minima 
become degenerate and sX T < T^, {rjt^cT{\ turns out to 
be the true vacuum (c/. Fig. [2]). At some temperature 
To ^ Tc the false vacuum decays into the true vacuum 
via bubble nucleation, described in section IV Bl Fig. [3] 
shows this evolution along the straight-line path connect- 
ing the two minima, as universe cools down. In Fig. ([4]), 
we show the evolution of the order parameter, the SM 
Higgs VEV r]. 




100 200 300 400 500 

(p (GeV) 



FIG. 3: The potential configuration along the straight-line 
path connecting the two minima, at various temperatures. 
At very high T the potential possesses only one minimum. 
At a lower temperature Ti a local minimum starts forming. 
At Tc < Ti the two minima become degenerate. At To < Tc 
tunneling to the true vacuum occurs. At T = the universe 
has settled in the true vacuum. The curves correspond to 
parameter set A of table HI 



B. Phase Transition through Bubble Nucleation 

The tunneling from the false to the true vacuum was 
calculated numerically using an approximation in which 
the bounce [s^ was assumed to lie along a straight line 
in the 2-d field space (?7, cr). Let (p be the field con- 
figuration along this path. At finite temperature, one 
looks for solutions of the Euclidean equations periodic in 
the "time" direction with period T~-'^[3^. In the high- 
temperature limit, the solution should be translationally 
invariant along the "time" axis, thus the dependence of 
on temperature disappears. The 0(3)-symmetric (in the 
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50 



FIG. 4: The evolution of the SM Higgs VEV rj with tempera- 
ture. At very high T, symmetry is restored: rj = 0. At a lower 
temperature a second order phase transition to ^ takes 
place, while still remaining in the false vacuum (dashed line) . 
At T = To, a first order phase transition brings the universe 
to the true vacuum (solid line). At T = 0, 77 = 246 GeV. The 
data corresponds to parameter set A of table H] 



spatial coordinates) solution will now obey the equations: 

(f(t) 2d(j) _dV d(t) 
dr^ r dr d6^ dr 



0, (j){oo) = 0, (33) 



where (j) — is the false vacuum. The decay rate per unit 
volume is 



T 



4 I S3(J)_Y^^ S3/T 

2ttT J 



(34) 



where we neglect the prefactor due to the change in the 
symmetry group (40| . Here Sslcj)] is the 3-dimensional 
action: 



(35) 



For the solution of eq. p3p we adopt numerical meth- 
ods, rather than using the well-known approximation 
schemes at the thin and thick wall limits [3^. The O4 
symmetric case at T = has been solved numerically 
in [Z^l. We do the same for the O3 symmetric equation 
and present here the results, used for the numerical 
estimates of table ID 

A potential of the form: 



1 



1 



1 



F(0) = -Af(T)202 _ _J(r)<^3 ^ _^(j.)^4 (3g) 

with > to ensure at least metastability at = 
0, encompasses all of the renormalizable potentials. For 
such a potential the transition between different regions 
depends on a single dimensionless parameter: 



(52 



(37) 



For tunneling to occur we must have k ^ while k ^ 
is required for the potential to be bounded from below. 



At finite temperature, the O3 symmetric action (|35|) is 
found to be: 



9M^ 



'3, thick 



125.8k- 
33k 



239.5k' 



1 -4k 



(38) 



with 53,thick ^ 13.72(1 -4k)2 and 53,thi„ = '^3^/" ^^X)- 
being the limits for the thick and thin wall approxima- 
tion, respectively [39j . 

The time needed for the universe to tunnel from the 
false to the true vacuum is estimated by setting T-tj^ « 1 
where tn ~ ^^"^ Hubble time in the radiation- 

dominated universe. Taking into account that during the 
electroweak phase transition T « 100 GeV, this yields: 
^ w 160, which defines the tunneling temperature Tq in 
the estimations of table |I1 If ^ 160, the tunneling 
occurs very quickly, when the universe is still hot. If 
160, the tunneling rate is too low and the universe 
remains at the false vacuum. 

When this condition is met, the universe tunnels from 
the false vacuum 0t to the true vacuum (^t. The en- 
ergy gained from the transition to a deeper minimum 
reheats the universe from the tunneling temperature To 
to a higher temperature T^. Since the expansion of the 
universe is much slower than the tunneling, the reheat- 
ing temperature is found by taking the energy den- 
sity to be constant p((/)f (To), To) = p{(j)t{Tr),Tr), where 
p(0,T) = f{(f>,T)+T s{^,T), with /(0,T) = Kff(0,T) 
the free energy density and s — ~df/dT the entropy 
density. In figH] the free energy density in the true and 
false vacuum are presented vs temperature. 



2x10' 



1x10' 



o 




-1x10" 



140 120 100 80 60 40 20 
T(Gey) 

FIG. 5: The free energy density vs temperature, at the true 
vacuum (lower curve) and at the false vacuum (upper curve). 
Significant reheating could occur if the universe cooled down 
to a low temperature before the tunneling took place. At 
the temperature at which the tunneling actually occurs, T ~ 
100 GeV, the reheating is not significant. 

Representative numerical estimations done using the 
above are presented in table HI The independent param- 
eters of the model were chosen to be A// , As, Ajyg, a, w 
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parameter sets 


A 


B 


Xh 


0.5 


0.6 


As 


0.6 


0.4 


Xhs 


0.025 


-0.02 


a 


2 


- 25 


(jj 


25 


90 




-200 


-220 




266 


220 


To 


168 


179 


Tr 


170 


182 




1.4 


1.5 


m*i,2(0) 


227, 247 


188, 302 




191, 139 


118, 101 




215, 201 


251, 143 


{P)/T 


0.84 


0.85 


f 


4.2 • 10"*^ 


4 • 10"® 


nis (keV) 


8.42 


8.85 



TABLE I: Representative parameter sets (see text for dis- 
cussion). The unit for all dimensionfull parameters is GeV, 
except for m^, which is given in keV. 



and the VEVs of the two Higgs bosons at zero tempera- 
ture r]o,ao, of which 7]q is kept fixed at 246 GeV. 

Tc, To, Tr stand for the critical, tunnehng and reheat- 
ing temperature respectively. Vc is the distance between 
the two vacua at the critical temperature and |f > 1 is 
the criterion for a strong 1st order phase transition. The 
parameter space of the potential ([4]) can provide for a 
variety of phase transition scenarios (1st order only, 2nd 
order only, 2nd order followed by 1st order). The parame- 
ter sets presented here fulfill the requirement for a strong 
1st order P.T. In A, a second order phase transition to 
non-zero VEV of the SM Higgs precedes the first order 
one to the true vacuum (fig. |4]). In B, no second-order 
phase transition occurs. 

"^1 2(^)1 "^1 2(^) stand for the Higgs mass eigenval- 
ues at the true and false vacuum respectively, at tem- 
perature T. For the parameter sets of table HI both of 
the Higgs modes decay after the transition to the true 

vacuum, since Tr > '"'"23^''"' 1 which is the temperature 
at which the decay rate appears to be maximal. Con- 
versely, decay in the false vacuum would require tunnel- 
ing temperature small enough. To < , in order for 
the Higgs bosons to have time to decay before the phase 
transition, and also sufficiently heavy Higgs eigenstates 
in the true vacuum, m}{Tr) > 2.3 T^, so that decays after 
the phase transition are suppressed by the low number 
density of Higgs bosons. No parameter sets satisfying the 
above were found. 

The values of the sterile neutrino Yukawa coupling / 
to the Higgs singlet, presented in table HI are obtained 
by requiring that sterile neutrinos make up all the dark 
matter, where now the details of the two-component de- 
cay, the phase transition and the decoupling of degrees of 



freedom were taken into account. The numerical results 
are consistent with the estimate of eq. (HU. The sterile 
neutrino mass = / • (Tq is then set to be in the keV 
range. 

VI. STERILE NEUTRINO PRODUCTION 
FROM OUT-OF-EQUILIBRIUM DECAYS 

Finally, we address the possibility of S decoupling early 
from equilibrium and decaying into sterile neutrinos out 
of equilibrium. This is the case if a, a; ~ and Xhs ~ 
10~^. Then, only a second order, rather than a first 
order, phase transition takes place and the S decays occur 
in the broken phase. 

The sterile neutrino population is again found from 
eq. (|10p . where now we need to first determine the out- 
of-equilibrium concentration of S bosons. 

The S boson number density Ns after decoupling, tak- 
ing into account the annihilations of S bosons to SM par- 
ticles, is given by [23j : 

NsjT) N^Tj) 1 
T3 Tf 2-rf/r ^ ^ 

If it were only for the SS XX annihilations, the S 
boson abundance would decrease at r ^ 00 to just half 
of its equilibrium value at freeze-out. However, the de- 
cay of S particles to sterile neutrinos will result in an 
exponential damping of the S boson abundance. In ad- 
dition, after H and 5* develop VEVs, S bosons will decay 
to SM fermions through the mixing with the SM Higgs. 
We can therefore ignore the 5* 5* — > XX annihilations and 
consider only the S ^ NaNa and S ff decays to de- 
termine 715 after freeze out. The kinetic equation for S 
bosons is: 

where 

^ \^ vn? j \m'Ax{\H,Xs)J 

takes into account the decay to all of the sterile neutrino 
species and SM fermions ^. Here A/ are the Yukawa 
couplings of the SM fermions to the Higgs doublet and 



We note in passing that sterile neutrinos with MeV< Ma < 
It. 2), produced in the S decays, can decay into three 
active neutrinos via mixing, in a tree-level process that involves Z 
exchange. Hence, we must take into account only the branching 
ratio of decay into the long-lived singlets with with keV-scale 
masses. 
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^f^^ , V is the mixing angle of the two Hiegs mass 

eigenstates, at the hmit Xhs '^h, A5 and a rj. Since 
S bosons hve in the electroweak scale, the main fermion 
decay mode will be the 6 b channel. For Xhs ^ 10~^, 
this dominates over the decays into keV sterile neutri- 
nos. As we will see below, Xhs ~ 10~® and Ma ~ keV is 
a self-consistent set of parameters for producing a suffi- 
cient amount of sterile neutrinos to make up dark matter, 
through out-of-equilibrium decays of S bosons. 
In terms of r = m/T and xs — Ps/T, one obtains: 



dns 
dr 



h^Mo 



=ns. 



(42) 



This yields: 

ns{xs,r) = 



1 f r + ^Jxl + r2 



where we set A = '}„ and we took ns to be the thermal 
equilibrium distribution function at T = Tf. 

Using ([42|) , eq. (fTO|l can be partially integrated to give 
the sterile neutrino distribution function, produced by S 
bosons decays after their freeze-out: 



0. . 
n [x,r) 



B 



I xsns{xs,rf)dxs 



4x 



xsnsixs,r)dxs 



where we set n 0(a;,rf) = and 



dr' 



(44) 



(45) 



is the branching ratio of 5* ^ NiNi to all other decays, 
with Ni being the lightest sterile neutrino. The last term 
in ((44|) vanishes at the limit r ^ 00, while the second 
term does not contribute to the total abundance, but 
only shifts the momentum distribution. The abundance 
of sterile neutrinos at any later time will be proportional 
to the amount of S bosons that have already decayed up 
to that time: 



Yf{r) = B[Ys{n)-Ys{r)] 



where 



Ys{r) 



45 



ns{xs.,r)x1dxs 



(46) 



(47) 



is the S boson abundance, with ns{xs,r) given by (|43p . 
The production rate dY^ /dr peaks at some rpi-od — 



m/Tpj-od- Since A determines how fast 5* bosons de- 
cay, Tprod depends on A but is effectively independent 
of r[. Given that Xhs ^ 10'^ is required for S bosons 
to be in equilibrium at early times, A receives a mini- 
mum contribution from the bb decay mode, through the 
mixing with the SM Higgs. Taking into account the b 
quark Yukawa coupling to the SM Higgs Af, ~ 2 • 10"^ 
and that Xh,Xs < 1, the decay into bb pairs ensures 
that A ^ 0.01, which results in 5* bosons decaying early 
enough, at rprod < 10, before the decoupling of the QCD 
degrees of freedom. Then, in eq. ([i7|) . 5, « 90 — 110, or 
^ « 25 — 33. The dependance of the final sterile neutrino 


abundance Y^ {00) /B on rf is shown in Fig|Sl 



(43) -i 




FIG. 6; The final sterile neutrino abundance {00) /B vs 
Tj. ~ 1 — 2 corresponds to out-of-equilibrium decays of 
S bosons, producing sufficient amount of dark matter. For 
> 3, the amount produced is insignificant in comparison to 
that produced during the in-equilibrium decays, occurring at 
r < 3, and can be ignored. 

The amount of dark matter produced from the out-of- 
equilibrium decay of 5" bosons is: 




(48) 



Early decoupling of S bosons 1 ^ ^-f ^ 3 implies 
Xhs — 10~^ (c/. Fig. [T|). The dominant decay mode is 
then bb pairs and the branching ratio of S decays into an 
nis ~ 3keV sterile neutrino, i.e. with / = (1 — 5) x 10^^, 
is S ~ 0.1 — 0.01. The amount of sterile neutrinos pro- 
duced by the out-of-equilibrium decays is then sufficient 
to constitute dark matter. 

The average momentum at r ^ 00 of the ster- 
ile neutrino population produced through the out-of- 
equilibrium decays is 



(p) ^ /r7 ^^'^^ /o°° dxx'^ns {x, r) 



t 

dxx'^ns{x,rf) 



(49) 



The variation of {(p)/T) with A is shown in Fig. [3 for 
various values of rf, where the redshifting factor ^^/^ 
(eq. (fT5|) ) has not yet been included. Then, for the set of 
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0.01 0.1 1 10 



A 

FIG. 7: Final average momentum (p)/T vs A cx /i^, for var- 
ious values of r;, stated next to each curve, (p) /T is further 
redshifted by ^^''^ ~ 3.2 with respect to the values presented 
in this graph. Small A implies delayed decays, resulting in 
warmer dark matter. In the out-of-equilibrium decay sce- 
nario, A ~ 0.1 results in production of sufficient amount of 
dark matter. 



three sterile neutrinos below the electroweak scale, and if 
the two heavier ones are nearly degenerate in mass [s^]. 
This scenario is different from the more commonly dis- 
cussed thermal leptogenesis in that neutrino oscillations, 
not decays, are responsible for the change in the lep- 
ton number of the plasma. Active neutrinos (in equi- 
librium) can oscillate into the sterile neutrinos (out of 
equilibrium), and CP violation in the neutrino mass ma- 
trix could make the the net lepton number of the out-of- 
equilibrium sterile neutrinos non-zero. The excess lepton 
number remaining in plasma is partially converted into 
the baryon number by sphalerons [sO] . 

An alternative possibility exists if the phase transition 
is strongly first-order, which is quite likely in the model 
with the singlet, as discussed in Refs. [35', 36, In this 
case the standard electroweak baryogenesis j44|i] can take 
place in the course of this first-order phase transition. 
The model of eq. ^ can easily be modified to include a 
sufficient amount of CP violation: all that is required for 
a successful baryogenesis is to include the second Higgs 
doublet [H]. 



parameters discussed above A w 0.01 — 0.1 and the pro- 
duction peaks at rprod — 3. Thus, the average momen- 
tum at the current temperature, for decoupling around 
rf ~ 1 — 2, can be as low as (c/. Fig. [T]): 



M 

T 



2.5 



= 0.8. 



(50) 



Finally, we return to the in-equilibrium decay scenario. 
If rf > 3, 5 bosons decay primarily while still in equilib- 
rium. However, the amount of dark matter produced 
will be supplemented by any additional sterile neutrinos 
produced after S bosons come out of equilibrium. This 
again is given by eq. (|48|) . where now both Yg/B and B 
take lower values because of the increase in r{ and Xhs 
{cf. Fig. [6]). Thus, the amount of dark matter produced 
is determined only by a fraction of what was produced 
before the freeze-out, and our results from the previous 
section remain valid. 



VII. BARYOGENESIS 

The model under consideration and its minimal mod- 
ifications offer at least two scenarios for generating the 
baryon asymmetry of the universe below the electroweak 
scale. One possibility is that the baryon asymmetry could 
arise from the low-scale leptogenesis if there are at least 



VIII. CONCLUSION 

The inclusion of singlet fermions (right-handed neutri- 
nos) in the Standard Model is the usual way to gener- 
ate the observed neutrino masses [ij. In contrast with 
many other models, we assume that both Dirac and Ma- 
jorana neutrino masses are generated via the Higgs mech- 
anism [20l| . The immediate advantage of this model is the 
possibility to produce dark matter, in the form of sterile 
neutrinos, which is cold enough to satisfy the bounds on 
the small-scale structure and the bounds from X-ray ob- 
servations, while explaining the pulsar kicks at the same 
time Q- The same sterile neutrinos can play an im- 
portant role in the formation of the first stars |jd| . We 
have considered different ways in which the dark-matter 
sterile neutrinos can be produced from the Higgs decays 
in the early universe. If the production from the Higgs 
decays dominates over the production by neutrino oscil- 
lations, the resulting dark matter population is colder 
than in the Dodclson-Widrow case for the same mass. 
The Higgs structure of the model has important implica- 
tions for the collider physics an can be probed at the at 
the Large Hadron Collider and a Linear Collider [s^, ■ 
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